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ABSTRACT 

Outflows form an integral component in regulating the gas cycling in and out of galaxies, 
although their impact on the galaxy hosts is still poorly understood. Here we present an anal- 
ysis of 405 high mass (log M./Ms > 10), star-forming galaxies (excluding AGN) with low 
inclinations at z ~0, using stacking techniques of the Nap 445889,5895 À neutral gas tracer 
in IFU observations from the MaNGA DRI15 survey. We detect outflows in the central re- 
gions of 78/405 galaxies and determine their extent and power through the construction of 
stacked annuli. We find outflows are most powerful in central regions and extend out to ~1Re, 
with declining mass outflow rates and loading factors as a function of radius. The stacking of 
spaxels over key galaxy quantities reveals outflow detections in regions of high Xsrg (20.01 
Moyr-!kpc^?) and Ey, (2107 Mokpe~’) along the resolved main sequence. Clear correla- 
tions with Xsrg suggest it is the main regulator of outflows, with a critical threshold of ~0.01 
Moyr^!kpc ^? needed to escape the weight of the disk and launch them. Furthermore, mea- 
surements of the Hô and D,,4000 indices reveal virtually identical star formation histories 
between galaxies with outflows and those without. Finally, through stacking of H 121 cm ob- 
servations for a subset of our sample, we find outflow galaxies show reduced H : gas fractions 
at central velocities compared to their non-detection control counterparts, suggestive of some 
removal of Hı gas, likely in the central regions of the galaxies, but not enough to completely 
quench the host. 
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1 INTRODUCTION 


The evolution of galaxies is dictated primarily by the availability 
and regulation of cold gas available for star formation, both at low 
and high redshift (Saintonge et al. 2013). In this scenario, a galaxy 
fills its gas reservoir through accretion, converts the gas into stars 
over a free fall time, and ejects the subsequently metal-enriched 
gas out of the galaxy disk (Davé et al. 2012; Lilly et al. 2013). 
Such a picture is consistent with recent molecular and atomic gas 
observations of representative galaxies at z ~0, which show sig- 
nificantly elevated gas fractions and star formation efficiencies for 
star-forming galaxies along and above the so-called “star-forming 
main sequence" compared to their passive counterparts (Saintonge 
et al. 2017; Catinella et al. 2018). 

However, the conversion of gas into stars is an inefficient pro- 
cess, with scaling relations showing that only ~1-10% of the gas 
gets converted into stars at a given time (Kennicutt 1998; Usero 
et al. 2015; Bigiel et al. 2016). Galactic-scale outflows, launched 
through the combined energy and momentum delivered by intense 


* E-mail: guidorb@astro.ucla.edu 


© 2019 The Authors 


star formation or a supermassive black hole (an active galactic nu- 
cleus; AGN) are thought to contribute to this inefficiency through 
their potential to regulate and/or halt the star formation in the disks 
of galaxies, thereby providing the feedback necessary to deliver 
galaxies from the star-forming main sequence to the red cloud of 
passive galaxies. AGN feedback is thought to be the dominant regu- 
lating process for high mass galaxies at later times, whilst starburst- 
driven outflows are thought to dominate low mass galaxies at early 
times. 


However, outflows are still poorly understood. Recent efforts 
aimed at understanding them have made progress by placing con- 
straints on their prevalence, integrated properties, and the under- 
lying correlations with their galaxy hosts through observations of 
relatively large samples of galaxies (Rupke et al. 2005; Chen et al. 
2010; Sugahara et al. 2017; Roberts-Borsani & Saintonge 2018). 
Such studies have found that 100-1000 km s^! outflows are com- 
mon at all epochs among star-forming systems and AGN (Veilleux 
et al. 2005; Rupke et al. 2005; Feruglio et al. 2010; Chen et al. 
2010; Coil et al. 2011; Cicone et al. 2016; Cazzoli et al. 2016; 
Roberts-Borsani & Saintonge 2018) and particularly at high mass 
(log M./Mo > 10) for normal galaxies of the local Universe along 
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the main sequence when traced in absorption (Chen et al. 2010; 
Roberts-Borsani & Saintonge 2018; Concas et al. 2019; Sugahara 
et al. 2017). Estimates of a mass outflow rate in star-forming galax- 
ies, whilst extremely challenging and uncertain, have revealed val- 
ues similar to the SFR of the galaxy, suggesting that more intense 
star formation expels higher quantities of gas (Rupke et al. 2005; 
Rubin et al. 2014; Cazzoli et al. 2016; Roberts-Borsani & Sain- 
tonge 2018). Similar correlations are found between the mass out- 
flow rates and AGN luminosities for galaxies with an active nu- 
cleus, suggesting more extreme AGN expel gas at significantly 
higher rates (Sturm et al. 2011; Veilleux et al. 2013; Cicone et al. 
2014; González-Alfonso et al. 2017; Fluetsch et al. 2019). 

Such constraints have proved extremely valuable, however 
they are generally based on spectroscopic observations with limited 
and incomplete spatial coverage (e.g., long-slit or 3" fiber) and as 
such may not be representative of entire galaxy systems. The need 
for more detailed and representative galaxy observations has seen 
the advent of large integral field unit (IFU) spectroscopic surveys 
such as CALIFA (Sánchez et al. 2012), SAMI (Croom et al. 2012) 
and MaNGA (Bundy et al. 2015), which have heralded in a new era 
for statistical studies of representative galaxies and their outflows 
at parsec (pc) and kiloparsec (kpc) resolution. The arrival of such 
large IFU surveys, in conjunction with powerful instruments such 
as VLT/MUSE, have added valuable constraints on the prevalence 
of ionised-gas outflows in normal, individual galaxies, their driv- 
ing properties, star-formation histories and AGN contributions on 
a per-spaxel basis (e.g., Ho et al. 2016; López-Cobá et al. 2019; 
Rodríguez del Pino et al. 2019). Despite such progress, several 
outstanding questions remain: whether outflows can quench their 
hosts, either entirely or in specific regions, remains an open ques- 
tion and the primary driving properties in star-forming galaxies has 
yet to be fully determined. 

In this paper, we expand on recent IFU work by using one 
of the largest IFU samples to date of star-forming galaxies ob- 
served as part of the MaNGA survey containing signatures of neu- 
tral gas (Nap) outflows, with the goal of determining through stack- 
ing techniques (1) the radial extent of outflows and their power, (ii) 
the main driving mechanism of outflows and the Nap absorption 
doublet, and (iii) the average star-formation histories and H1 gas 
fractions of outflow galaxies compared to those without. The paper 
is organised as follows. In 82 and $3 we present our sample selec- 
tion, relevant data sets and stacking techniques using the MaNGA 
DR15 survey, respectively. $4 presents an analysis of the radial ex- 
tent of outflows and their power, whilst $5 and $6 explore the pri- 
mary driving properties of outflows. We present the average H1 gas 
fractions of outflow galaxies and a control sample of galaxies with- 
out outflows in §7. Finally, we provide a discussion and the star 
formation histories of outflow galaxies in §8 and present our con- 
clusions in $9. Throughout the paper we adopt a ACDM cosmol- 
ogy with H9270 km/s/Mpc, Qm 20.3, and O4 =0.7, and assume a 
Chabrier IMF. 


2 DATA AND SAMPLE 
2.4 The SDSS-IV/MaNGA Survey 


The SDSS-IV/Mapping Nearby Galaxies at APO (MaNGA; Bundy 
et al. 2015) survey is an IFU survey of the local Universe with the 
aim of obtaining spatially-resolved spectra for an unbiased sam- 
ple of ~10,000 galaxies from the NASA-Sloan Atlas (NSA) cata- 
log, with log M./Mo 29 at 0.01 <z<0.15 by 2020. Situated on the 


SDSS 2.5m telescope at Apache Point Observatory, the MaNGA 
instrument offers 29 differently-sized IFUs, each of which consists 
of a set of optical fibers grouped to form a hexagon and fed into the 
two BOSS spectrographs (Smee et al. 2013). These fiber bundles 
are split into 17 sets which can be used to observe chosen targets 
at a given time, with the remaining 12 used for flux calibration and 
additional single fibers for sky subtraction (Drory et al. 2015). The 
bundles cover up to 1.5 Re and 2.5 Re for the targeted galaxy sam- 
ple and the resulting spectra have a wavelength coverage of 3600- 
10,400 À with spectral resolution R~2000 (Wake et al. 2017). On 
10th December 2018, the MaNGA survey released its 15th Data 
Release (DR15; Aguado et al. 2019), along with data products such 
as maps of derived galaxy properties and continuum fits from the 
Data Reduction Pipeline (DRP; Law et al. 2016) and Data Anal- 
ysis Pipeline (DAP; Westfall et al. 2019), as well as the Marvin 
Python package (Cherinka et al. 2019) to view and analyse them. 
Additionally, the release also comes with new value added catalogs 
of spatial and integrated properties such as the Pipe3D (Sánchez et 
al. 2018) and H1 follow up observations (Masters et al. 2019). In 
this study we make use of this release, which comprises completed 
observations for 4,656 galaxies. 


2.2 Sample Selection 


Motivated by recent results from Roberts-Borsani & Saintonge 
(2018), we begin our selection by performing cuts on galaxy stel- 
lar mass (log M./Mo >10), SFR (log SFR/Mo yr! »-2.332 (log 
SFR/Mo yr!) + 04156(log SFR/Mo yr7!)? - 0.01828 (log 
SFR/Mo yr!» - 0.4, corresponding to a rough lower limit of the 
galaxy main-sequence defined by Saintonge et al. 2016) and incli- 
nation (i <50°, derived from the galaxy's r-band axis ratio, b/a) - 
three key galaxy properties known to influence the detection rates 
of neutral gas outflows. The stellar masses and SFRs used for this 
selection are taken from the Pipe3D catalog, which for each galaxy 
derives a SFR based on the integrated Ha luminosity, and the axis 
ratio of the galaxy is taken from the NASA-Sloan Atlas catalog. Fi- 
nally, given that broad line regions in AGN can cause overestima- 
tion of Ha-derived SFRs and/or mimick the presence of outflows 
in ionised gas tracers, we choose to identify and remove objects 
with AGN signatures in their central regions using a Kauffmann 
et al. (2003) BPT diagnostic. Motivation for this also comes from 
several recent studies (e.g., Sarzi et al. 2016; Roberts-Borsani & 
Saintonge 2018; Concas et al. 2019) which have demonstrated the 
limited influence of a weak AGN in the normal galaxy populations. 
Our resulting sample consists of 422 star-forming galaxies. To en- 
sure our sample is not contaminated by the presence of mergers, 
stars, or pointing offsets, we visually inspect the MaNGA footprint 
image of each galaxy: 17 galaxies galaxies fail these criteria and 
are removed from the sample. Our final sample, therefore, consists 
of 405 star-forming galaxies, whose position on the SFR-M. plane 
is shown in Figure | along with histograms of their SFRsyq and 
stellar masses. The sample spans virtually the full MaNGA red- 
shift range, with a median redshift of 720.04 which corresponds to 
a pixel sampling of ~0.4 kpc and an effective spatial resolution of 
FWHM-2 kpc. 
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Figure 1. The SFR-M. plane and density contours of the full MaNGA 
DRIS5 sample and selected sample for this study. Gray dots represent the 
405 galaxies found using our selection criteria (blue dashed lines), whilst 
the blue points mark galaxies found to have outflows in their central regions 
and magenta points galaxies found to have inflows. The orange solid and 
dashed lines mark the star-forming main sequence defined by Saintonge et 
al. (2016) and its lower and upper limit. Histograms of the selected galax- 
ies’ SFRs and stellar masses are shown to the top and right sides of the main 
plot. 


3 STACKING PROCEDURES AND ANALYSIS 
3.1 Maps of Galaxy Properties 


As a first step in our analysis, we create maps of spatially-resolved 
galaxy properties (i.e., SFR, sgg, Mx, X., Ay and D(4000)) using 
the MaNGA DR15 Pipe3D IFU maps for each galaxy in our se- 
lected sample. We begin by using the spatially resolved Ha and HG 
emission, for which the underlying continuum has been removed 
by the MaNGA DAP, in order to derive a Balmer decrement for 
each spaxel, which we translate into an Ay and Agr, magnitude, 
assuming an intrinstic ratio of Ho/H22.68. The maps of Ha are 
subsequently corrected for dust and converted to a luminosity using 
luminosity distances derived with an Ha redshift and the assumed 
cosmology. The Hæ luminosities are subsequently converted to a 
Xsrpg using a Kennicutt (1998) prescription converted to a Chabrier 
IMF (SFR [Mo yt! ] Lg /[2.1x 107! erg/s]) and the physical area 
probed by each 0.5" spaxel. Not all spaxels, however, are appropri- 
ate for analysis and we therefore ensure a quality control by apply- 
ing the MaNGA bitmask flags and require the following criteria for 
science use: 


e A line S/N>3 for Ho, Hf, [On]45006, [N1]46583 and 
[O 146300. 

e A BPT “star-forming” nature determined by the combination 
of the [N 1]46583 Kauffmann et al. (2003) and [O 1]46300 Kewley 
et al. (2006) prescriptions. 

e Anr-band S/N>5 to guard against spaxels with very little con- 
tinuum signal. 


The redshifts for each spaxel are derived from the Talbot et al. 
(2018) value added catalog!. 15 galaxies do not have determined 


l https://www.sdss.org/dri5/data 
access/value-added-catalogs/?vac id- 
manga-spectroscopic-redshifts 
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spaxel redshifts and we exclude these from our analysis, resulting 
in a sample of 390 galaxies. In total, this results in 276,619 science 
spaxels and an example of the MaNGA maps for a representative 
galaxy in our sample is shown in Figure 2. We find that the planes 
of galaxy properties as traced by the full spaxel sample are fully 
sampled by each of the galaxy in our sample, ensuring stacks over 
any particular region of parameter space include virtually the full 
sample of galaxies. 


3.2 Stacking Procedure and Outflow Modelling 


Throughout this study we make use of stacking approaches to con- 
struct high S/N composite spectra with which to analyse outflow 
properties as a function of galactocentric radius and a variety of 
galaxy or spaxel properties. The spaxels are first divided into a set 
of bins according to a given property and each spectrum in an asso- 
ciated bin is subsequently corrected for foreground galactic extinc- 
tion using the associated Schlegel et al. (1998) E(B-V) values and 
an O'Donnell (1994) Milky Way extinction curve, after converting 
the wavelength arrays of the spectra to air wavelengths. Each spec- 
trum is then shifted to the rest-frame, before being interpolated over 
a common wavelength grid. 

We subsequently mask all flux points in a given spectrum 
that are deemed unfit for science by its associated mask array and 
normalise the spectrum by the median flux between 5450 À and 
5550 À (since this region is free of absorption and emission lines), 
thereby giving equal weight to each spectrum. The final spectrum 
is then taken as the mean over all the normalised spectra and the 
associated uncertainties are derived by adding in quadrature the 
bootstrapped sampling uncertainties and the mean flux uncertain- 
ties of each individual spectrum within the stack. In order to model 
outflow quantities associated with Nap, we use the same approach 
as described in Roberts-Borsani & Saintonge (2018) to model the 
line and use the blueshifted absorption component to derive outflow 
quantities. We refer the reader to their study for details, however we 
summarise here the basic procedure of the approach: each Na p pro- 
file is first fit with an analytical model (where the intensity of Nap 
is based on a velocity-independent covering fraction, linewidth, 
column density and a possible velocity offset. See Roberts-Borsani 
& Saintonge 2018 for details) of a single fixed and single offset 
profile to determine whether a flow is present, and subsequently 
characterised with a three-component profile (consisting of a sys- 
temic component, blueshifted absorption and redshifted emission) 
if an outflow is deemed present. We use the same priors discussed in 
Roberts-Borsani & Saintonge (2018), however given the low cov- 
ering fractions found in recent studies we limit the allowed values 
to |Cr| «0.5. 

The mass outflow rate can be estimated through knowledge of 
the local covering fraction (or clumpiness of the gas), Cf, hydrogen 
column density, N(H), wind velocity, Vout, and an assumed geom- 
etry and radius. Whilst most absorption studies of outflows derive 
mass outflow rates based on a spherically symmetric thin shell ge- 
ometry emanating from the centre of the galaxy, the resolved nature 
of the data sets use here means our stacks do not necessarily follow 
a well defined geometry. This is particularly relevant when stacking 
spaxels from different regions of a galaxy. Whilst one could derive 
rates based on the area probed instead, such an approach would 
result in direct correlation between the mass outflow rate and the 
number of spaxels in each stack, thereby removing much of any 
correlation with galaxy property. As such, we opt only to assume 
a radius of the outflowing gas of 1 kpc - where the outflowing gas 
is unlikely to be collimated by the disk - and not to assume a ge- 
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Figure 2. An example of IFU MaNGA galaxy property maps using the Marvin Python package. From left to right, top to bottom: The SDSS image of an 
example galaxy with the MaNGA footprint, the velocity of Ha relative to the systemic, the flux of Ha, the flux of H£, the Balmer decrement, the D(4000) 
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ometry. The advantage of this is that our absolute values are less 
subject to uncertain geometrical assumptions and are almost en- 
tirely derived from our fitted parameters, whilst the disadvantage 
is that the true values are likely to be higher than what we present 
here. Thus, our expression for a mass outflow rate becomes: 


Mout = Cr umy N(H) v r, (1) 


where each term is as defined above, and my is the mass of hy- 
drogen, with an assumed correction for the abundance of Helium, 
y=1.4. Unless explicitly stated otherwise, we assume this general 
formulation for the mass outflow rate throughout all sections of 
this paper, and the assumption of r=1 kpc also remain unchanged 
throughout. 


4 THE GALACTOCENTRIC PROFILE OF OUTFLOWS 


We begin by stacking spaxels as a function of deprojected galac- 
tocentric radius (i.e., accounting for the galaxy's inclination and 
rotation on the plane of the sky) for our sample of 390 galaxies, in 
order to create multiple independent annuli for each of our galaxies 
which probe the full range of the galaxy disk. The first annulus is 
centred at the centre of the galaxy disk and subsequent annuli ex- 
tend outward from the border of the inner annulus. Each annulus 
contains spaxels within a full width of 0.25 Re from its centre, and 
the distribution of annuli probe from the centre of the galaxy out to 
^2 Re. For each galaxy, all relevant spaxels are deprojected from 
the position of the central spaxel using the galaxy's position angle 
(PA) and inclination: 


RA geproj = (RA - RAcentzal) © cos(PA) + (DEC — DEC central) - sin(PA) (2) 


—(RA - RA central) + Sin(PA) + (DEC — DEC entra) : cos(PA) 
cos(i) 


DECateproj = 


(3) 


The deprojected radius between a given spaxel and the central 


spaxel is then taken to be rgeproj = [RAR oj + DEC $ eproj’ and the 
mean spectrum and quantities quoted here are those over the cen- 
tral spaxels. The motivation here is to determine and select galaxies 
that display signatures of outflows, in order for them to be used in 
subsequent analyses where we can tie global galaxy properties to 
outflow properties. Since the signature of blueshifted absorption 
can easily be “diluted” by strong systemic absorption in a stack, 
we wish to maximise our chances of outflow detection and char- 
acterisation by preselecting galaxies with outflows in their central 
regions. Thus, to determine whether a gas flow is present in the cen- 
tral region (R<0.25 Re), we compare Nap fits of a fixed systemic 
component to a single blueshifted or redshifted absorption com- 
ponent. For this we use a Bayesian Information Criterion (BIC), 
which penalises for extra free parameters and favours lower val- 
ues. We derive ratios (K=BIC¢xeq/BICgoy,) of the two models to 
account for the extra free parameter (i.e., a velocity offset) in the 
shifted model and determine a flow detection if the ratio is greater 
than unity and a minimum blueshift velocity is found. Furthermore, 
several recent studies have shown that Sodium excess can some- 
times be seen in emission (Chen et al. 2010; Roberts-Borsani & 
Saintonge 2018; Concas et al. 2019), with the rate of outflow de- 
tections decreasing rapdily as the net profile changes towards emis- 
sion. Given that the nature of the Nap emission is still poorly under- 
stood, we opt to analyse only galaxies which show net absorption at 
their centres, and require an absorption depth-to-noise (D/N) ratio 
>3 in the central region, with a minimum blue shifted of 15 km s^! 
if D/N»5 and 20 kms?! if 3<D/N<5. In total, 92 galaxies satisfy 
these criteria, with 78 objects displaying outflow detections and 14 
displaying inflow detections. The remaining 298 galaxies either do 
not show sufficient absorption (^ 80406), with a typical (median) D/N 
ratio of 1.62, do not show blueshifted absorption (~18%), or have 
insufficient spatial resolution to stack within 0.25 Re (~2%). We 
report the flow galaxies and their main properties in Table A1 and 
illustrate the typical D/N ratios in Figure 3. 
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Figure 3. The typical depth-to-noise ratios of Na D in the central («0.25 Re) 
regions of the sample of MaNGA galaxies used in this study. The gray his- 
togram represents non-detections, whilst the blue histogram and the ma- 
genta histgram represent outflow and inflow detections, respectively. Pos- 
itive values reflect profiles in absorption, whilst negative values represent 
profiles in emission. The dashed black line marks our minimum D/N thresh- 
old of >3 to determine detections of inflows and outflows. 


4.1 Average Central Profiles 


We begin by examining and comparing the mean spectra of galax- 
ies with and without flow detections in the inner 0.25 Re of their 
effective radii. We adopt a Monte Carlo approach to constructing 
the mean central stacks over samples of outflow and inflow galax- 
ies, using a control sample of non-detection galaxies for compari- 
son: a random sample of 10 outflow (or inflow) galaxies is selected 
and for each galaxy we randomly select a non-detection *coun- 
terpart", defined as such if it lies within +0.2 dex in stellar mass 
and SFR, before averaging the central spectra. This is repeated 
100 times and the final spectrum is simply the mean over all it- 
erations, with the stack flux errors derived from a combination of 
individual flux uncertainties and bootstrapping errors. The results 
of these are shown in Figure 4, where we show the differences in 
spectra between galaxies with outflows, inflows, and those with- 
out, along with some key galaxy properties. From the left panel, 
we note that the optical spectra over ~3700< 2 «8000 À for out- 
flow and non-detection galaxies are virtually identical, with minor 
differences only in the strength of the emission lines and depth of 
the absorption lines. The similarity is less pronounced for inflow 
galaxies, however, which display a considerably redder continuum 
slope than their non-detection counterparts. To quantify these dif- 
ferences, we take the mean continuum fluxes between 4500 < A < 
4800 A and 5500 < a « 5800 A - roughly equivalent to the B 
and V bands - and find B — V colours of 0.015+0.039 mag and 
0.003+0.041 mag for outflow galaxies and their associated non- 
detections, respectively, and 0.132+0.045 mag and 0.001+0.041 
mag for inflow galaxies and their associated non-detections, re- 
spectively. In the case of outflow galaxies and both non-detection 
spectra, the spectra display virtually no B — V excess, suggestive 
of a flat spectrum. However, this contrasts with the spectra of in- 
flow galaxies which display a more significant B — V excess (~0.12 
mag redder than the outflow galaxies), and fitting straight lines to 
the pseudo-photometry reveals an inflow galaxy spectral slope ~9x 
greater than that of the outflow spectra. 

A much starker contrast, however, appears when comparing 
the average Nap profiles (middle panel) profiles. For Nap we note 
a large difference in the total equivalent width (EW) of the pro- 
file, with ~0.5A for outflow detections, ~0.95 for inflow detec- 
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tions, and -0.1 A for non-detections. A visual comparison of the 
profiles shows that both inflow and outflow detections are charac- 
terised by significant and unambiguous absorption. The Na 145889 
line has an intrinsic line depth twice the size of its Na L45895 red- 
shifted counterpart (Morton 1991), however the flow profiles dis- 
play a clear asymmetry: the outflow spectrum is blueshifted with 
respect to the systemic wavelengths and the inflow spectrum, while 
not as significantly redshifted, shows nearly equal line depth for 
the two Nap lines which can arise due to an additional redshifted 
component. Both profiles are well-characterised by an additional 
blueshifted or redshifted component offset by |Av| «120 km s7! 
and |Av| ^85 kms! , respectively. The average non-detection spec- 
tra, on the other hand, are characterised by virtually no absorption 
and even display net emission. In stacking the full non-detection 
sample, we find this emission becomes even more pronounced. 

Although less pronounced, similar comparisons are found 
with the ionised gas traced by Ha-[N n] (right panel): the aver- 
age non-detection spectrum is characterised by a narrow profile 
that is less luminous than its detection counterpart, by a factor of 
0.7. In contrast to the outflow observed in the mean Nap pro- 
file, a significant broadening of the ionised gas is not immediately 
obvious, however a BIC ratio between a one- and two-component 
fit strongly favours the latter fit, with an outflow velocity of ~130 
km s7!, similar to that observed with the neutral gas and suggestive 
of the emission emanating both from the disk of the galaxy and a 
broader outflowing component. 

To gain some indication of the driver for the different pro- 
files, we compare the mean values of several galaxy properties 
likely to be important in determining both the shape of the Nap 
and Ha+[N n] profiles and flow presence. Nap requires significant 
amount of shielding from dust to survive, and dust obscuration has 
been found to correlate with the EW of the line (e.g., Veilleux et al. 
1995; Heckman et al. 2000; Chen et al. 2010; Rupke, & Veilleux 
2013; Roberts-Borsani & Saintonge 2018). Here we find average 
Ay magnitudes of ~1.4-1.9 mag (detections) and ~1.2 mag (non- 
detections), suggesting at first glance somewhat of a limited impact 
in regulating the shape of the two flow profiles. Such a conclu- 
sion may come as a surprise given previous results in the literature, 
however several reasons exist to explain this. The first is simply that 
dust is not the only regulating property of the EW of Nap, since in- 
clination and gas column densities will also impact the shape of the 
profile: indeed, Chen et al. (2010) find a strong correlation between 
NaD EW and inclination for SDSS galaxies which probe similar 
nuclear areas to our samples of MaNGA galaxies. With inclinations 
lower than i~40°, they observe systemic Nap transitioning from 
absorption into pure emission, possibly due to smaller column den- 
sities of gas probed along the line of sight through the thickness of 
the disk and difficulties of the continuum-fitting code to reproduce 
such low column densities. If we compare our results to theirs, we 
note that the median inclination of our high mass MaNGA galax- 
ies displaying net emission is significantly lower at ~35°, enough 
to push the Nap profile and EW into net emission. Finally, we do 
point out that even over the small range of mean Ay values, the 
relation between Na p EW and Ay is still present when comparing 
the values to the size of the Nap profiles in the middle panel of 
Figure 4. 

Similar observations are found for a dependence on the 
D(4000) break, which varies only by ~0.1 between the spectra, 
suggesting limited impact by the age of the underlying stellar pop- 
ulation. The difference in Xspr, on the other hand, we find to 
be much more significant: for non-detections, the average value 
(0.048+0.001 Moyr !kpc?) is lower compared to the average val- 
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ues of detections, which are both higher than >0.1 Moyr7!kpe~2 
(0.116+0.001 and 0.126+0.002 for outflows and inflows, respec- 
tively). Although the differences between the average values is 
not, at first glance, particularly large or significant, it is their ab- 
solute values which is of greater interest: the Usprs of the flow 
spectra are consistent with (and above) what is generally invoked 
as a minimum threshold for outflow activity seen in absorption 
(0.1 Moyr !kpc ?; Heckman 2002). Finally, we also find elevated 
stellar mass surface densities for the detection spectra, compared 
to the non-detections and this is most likely the result of a pos- 
itive correlation with Xspg (1.e., the “resolved” Main Sequence; 
e.g., Sánchez et al. 2013; Cano-Díaz et al. 2016; González Del- 
gado et al. 2016; Abdurro'uf, & Akiyama 2017; Hsieh et al. 2017; 
Maragkoudakis et al. 2017). We attempt to disentangle these intrin- 
sic correlations and explore further the main regulating properties 
of the outflows in Section 6. 


4.2 The Radial Extent of Outflows and Their Properties 


Outflows in star-forming galaxies are typically thought to display 
a degree of collimation in their geometry, most likely attributed to 
the underlying disk pressure gradient resulting from the interplay 
between gravity and the outward momentum from stellar feedback. 
Here, we test this picture by looking at the extent and properties of 
outflows over a range of galactocentric radii, R/Re. This is shown 
in Figure 5, where we plot the evolution of the ISM Nap EW (left 
panel), mass outflow rate (middle panel), and mass loading factor 
(right panel) as a function of radius. For reference, we also compare 
these to what would be measured by an SDSS 1.5’’-radius fiber. 

In the left panel of Figure 5 we plot the EW of the Nap resid- 
ual - measured directly from the spectrum itself - for stacks con- 
taining only the outflow galaxies (circles) and stacks with both the 
non-detection and outflows galaxies (diamonds). For the former, 
we observe a rapid decrease in EW as a function of radius, with a 
roughly exponentially declining profile which begins to flatten out 
slightly at 20.75 Re. The central regions of the galaxy («0.5 Re) 
are dominated by absorption, although this quickly transitions to 
net emission outward of —0.5 Re. A 1.5" stack of the detection 
galaxies shows that SDSS measurements would be consistent and 
similar to what we measure here. Although we cannot compare di- 
rectly, we note that this is similar to the distribution of molecular 
gas (i.e., cold gas) in disk galaxies traced by CO (Bigiel et al. 2008; 
Leroy et al. 2008; Schruba et al. 2011; Bigiel, & Blitz 2012), where 
observations of molecular gas as a function of galactocentric radius 
have found a roughly exponential decrease. Outflow detections are 
observed only out to ^1 Re, with the majority of the blueshifted 
absorption occurring within 0.5 Re and the outer regions detected 
primarily due to the imprint of redshifted emission that backscatters 
off of the receding side outflow. We note that some degree of bias 
may exist in these observations, since the detection of blueshifted 
absorption is contingent upon observing Nap in absorption, which 
in turn is partially regulated by the strength of the underlying con- 
tinuum. Since this latter consideration is driven by the surface den- 
sity of stars, which drops off steeply at higher radii, non-detections 
of outflows due to the lack of Nap absorption are likely to occur 
in regions of low Xy. As such, it is possible our detections are bi- 
ased towards the central regions of the galaxies where higher stel- 
lar mass surface densities are found, however our choice to stack in 
units of R/Re should provide some safety against such a bias. 

A milder trend is observed in stacks containing both the non- 
detection galaxies and outflow galaxies (diamond symbols in the 
left panel of Figure 5), albeit with some important differences: the 


strength of Nap emission increases slightly with radius, but in gen- 
eral maintains a much flatter slope compared to when only outflow 
galaxies are considered. Additionally, the added presence of the 
non-detection galaxies ensures the EW never reaches absorption at 
any radius and is observed only in emission; no detections of out- 
flows are observed in these stacks. The contrast between the two 
samples is particularly evident below ~1 Re, where outflow detec- 
tions are present, and becomes less great outward of this limit. Such 
a comparison highlights the ease with which signals of outflows in 
normal galaxies can be overlooked if not selected appropriately. It 
is unclear, however, what induces such a rapid change in the Nap 
profile. Although typically assumed to be an absorption transition, 
several studies using SDSS (Chen et al. 2010; Roberts-Borsani & 
Saintonge 2018; Concas et al. 2019) have also found Nap in emis- 
sion in stacked spectra. One obvious hypothesis is that SSP models 
at such low Nap EW are unable to accurately reproduce the stel- 
lar contributions and as such overestimate them, although there is 
a possibility that at outer radii there is insufficient gas and dust to 
produce significant Na p absorption and potential background Nap 
re-emission begins to dominate. 

We also wish to determine the radial evolution of the mass 
outflow rates in our stacks. To do this, we use Equation 1 to de- 
rive values of Mout for the detections presented in the left panel of 
Figure 5 and plot these as a function of radius in the middle panel. 
In these stacks, the adopted N(H) value is the mean value of the 
summed column density for a given annulus, over all galaxies in 
the stack. We immediately note an important decrease in Mout as a 
function of galaxy radius, similar to the trend observed for the Nap 
EW in the left panel. The outflow appears strongest in the central 
regions of the galaxy, with values starting at ~0.13 Mo yr! and de- 
creasing down to ~0.01 Moyr-!. A comparison with the 1.5" stack 
shows SDSS observations only probe the strongest parts of the out- 
flow, with a mass outflow rate of 0.14 Moyr7!, and likely overlook 
important contributions from the outer regions (i.e., out to ~1 Re). 
Specifically, an SDSS 3" fiber would miss ~50% of the total mass 
outflow rate integrated over all the MaNGA radial stacks, which 
we calculate to be Mout ~0.28 Moyr™!. The trend observed by the 
MaNGA points is well described by a linear fit, with a slope of ~-1: 


logio Mou, = (-1.18 + 0.12) - R/Re — (0.66 + 0.08) (4) 


Such a trend is perhaps not surprising: the outflows are se- 
lected to be star formation-driven and as such are likely to corre- 
late with star formation- and gas-dependent quantities, which are 
known to be mostly centrally concentrated (e.g., Bigiel et al. 2008; 
Ellison et al. 2018). 

The derivation of mass outflow rates in conjunction with the 
mean SFRs associated with each ring stack allow us to compute 
a resolved mass loading factor, 7, and to first order, determine the 
extent of the outflows’ potential for quenching. The mass loading 
factors are shown in the right panel of Figure 5. We observe a very 
similar trend to the mass outflow rate, with a rapidly decreasing 
logy as a function of radius, characterised by a slightly steeper 
slope of -1.4. 


logon = (-1.42 + 0.15) - R/Re — (0.22 + 0.10) (5) 


The values range from mass loading factors of 7 ~0.4 in the 
central regions of the outflow and decrease to 7 ~0.02 in the out- 
ermost regions, suggesting the potential for ejective quenching by 
the outflows is strongest in the central regions and not a galaxy- 
wide phenomenon - i.e., only in the central regions does the outflow 
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Figure 4. The average central 0.25 Re spectra of galaxies harbouring outflows (blue) or inflows (magenta), and the average spectra of an associated control 
sample of non-detection galaxies (light gray for inflows, dark gray for outflows). The 1o scatter is shown as the colour shaded regions around each spectrum. 
Left: the normalised optical spectra over virtually the full optical range. Middle: the continuum-normalised ISM residual of the Na p line for the four spectra 
in the left panel. Right: the continuum-normalised Ha+[N r1] emission for the four spectra in the left panel. The quantities and units quoted in the left panel are 
the stellar mass surface density (Mokpc ?), SFR surface density (Mo yr"! kpc?), dust (magnitudes) and 4000 A-break. The colours correspond to the different 
stacks (NB: given that the average spectra - and quantities - for the two non-detection control samples are essentially identical, we display only the values of 
the control sample associated with the outflow detections). 
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Figure 5. The evolution of the total Na D EW (left panel), mass outflow rate (middle panel) and mass loading factor (right panel) as a function of galactocentric 
radius. 1o error bars are included in each panel and in many cases are smaller than the symbols themselves. The Nap EW is separated out into stacks of 
spaxels from the outflow detection sample only (circles) and non-detection+outflow spaxels (diamonds). Filled symbols indicate the detection of outflowing 
neutral gas, whilst empty symbols indicate non-detections. For comparison, a 1.5” stack is also added to each of the plots, in order to gauge what would be 
measured by single-fiber SDSS surveys. In the middle and right plots, we add a best fit first order polynomial (dashed line) and its lo error (shaded region). To 
illustrate the impact of different geometric assumptions in the calculation of mass outflow rates and loading factors, in the right panel we plot the mass loading 
factor assuming no outflow geometry (dark blue points and orange SDSS point; as described in Section 3 and by Equation 1) and a spherically symmetric thin 
shell geometry of 4z (light blue points; these are simply the dark blue points multiplied by 47, see Section 4.2 for details). These serve as lower and upper 
limits of the mass loading factor, respectively. 


have any sort of potential to quench the galaxy host by removing bles" at larger radii (e.g., 5 kpc) above the disk. Thus, we multiply 


gas faster than the rate of star formation and halting star formation 
activity - although none of the values ever reach unity, suggesting 
quenching even in the central regions remains unlikely. Integrating 
over all outflow detections, we find the galaxy-wide mass loading 
factor is 7 «0.1, supporting this hypothesis. 


However, as mentioned in Section 3.2, our choice not to as- 
sume an outflow geometry for these outflow rates means we are, to 
a degree, likely underestimating the absolute values. As such, for 
comparison we calculate mass loading factors with outflow rates 
derived with an assumed geometry. We assume here the outflowing 
gas from each stack is coming from an individual “wind bubble", 
which we measure as a spherically symmetric shell subtending 47 
steradians, whose origin resides in the galaxy disk and extends out 
to 1 kpc. The base assumption here is that the large-scale outflow 
is formed via the superposition and collimation of such “wind bub- 
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the outflow rates derived with Equation 1 by 4z and show these as 
light blue points in Figure 5. The result shows mass loading fac- 
tors ~1.1 dex greater than our fiducial measurements, with a range 
0.2 7 £5.2, suggesting an enhanced possibility of first order ejec- 
tive quenching in the central regions of the galaxy and highlight- 
ing the importance of geometrical assumptions. We note that for a 
spherically symmetric thin shell geometry at 1 kpc, the two ranges 
of outflow rates presented here reflect the lower and upper limits of 
the model. 


It is unclear whether it is a radial dependence or the average 
galaxy properties which drive the main trends seen in Figure 5. In 
fact, by comparing galaxy properties associated with each radial bin 
to the evolution of the mass outflow rate, we find the latter quantity 
follows most closely the evolution of Xgpgg and £. This is shown 
in Figure 6, where we plot the normalised evolution of galaxy prop- 
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Figure 6. The normalised evolution of the geometry-independent mass out- 
flow rate compared to several key galaxy properties (i.e., Uspr, Xw,. Av. 
D(4000)), with error bars. Each quantity is normalised by the maximum 
value across galaxy radius and error bars of galaxy properties include - for 
each annulus - the standard deviation over all spaxels and all galaxies going 
into the stack. The radial evolution of the normalised mass outflow rate is 
most closely aligned with the normalised evolution of Xsrg and Xy,. 


erties and the mass outflow rate. However, from these stacks alone 
it is difficult to determine whether the selected properties or the ra- 
dial dependence are the primary regulators of the outflows, given 
that intrinsic correlations exist between galaxy radius and the cho- 
sen properties. We further inspect the correlation of (or lack of) the 
Nap EW and outflow properties with global galaxy properties in 
Section 6. 


5 THE RESOLVED gsrg-Xy, PLANE 


Since the discovery of the galaxy Main Sequence (MS) and the de- 
velopment of basic frameworks to describe a galaxy's position rela- 
tive to it, much work has gone into determining the prevalence and 
influence of outflows relative to the MS (e.g., Chen et al. 2010; Ru- 
bin et al. 2014; Cicone et al. 2016; Roberts-Borsani & Saintonge 
2018; Concas et al. 2019). The arrival of large IFU surveys such 
as MaNGA, SAMI and CALIFA, has also revealed a resolved MS 
(Sánchez et al. 2013; Cano-Diaz et al. 2016; González Delgado et 
al. 2016; Abdurro’uf, & Akiyama 2017; Hsieh et al. 2017; Maragk- 
oudakis et al. 2017), indicative of a link between small scale pro- 
cesses and the integrated properties of galaxies. As such, it is in- 
teresting to look at the prevalence and properties of outflows over 
such a resolved sequence. We therefore present a stacked analysis 
over the local MS in Figure 7, using spaxels from our sample of 
outflow galaxies. Here, the mean stacks and properties are taken 
over all spaxels in a bin, rather than over the galaxies going into 
the stack. From the left panel of Figure 7, we find a very similar 
trend of outflow prevalence and strength to what is found in in- 
tegrated analyses: outflows are found predominantly in regions of 
high Xspg and Xj, , with increasing strength (as traced by the out- 
flow EW and mass outflow rates) higher up the MS. Here we find 
that detections span a X. range of 7.5«log XJMokpc? «9.5 and 
Xgpg of 2«log Xgpgg/Moyr !kpc ? «0, in agreement with lim- 
its found by Roberts-Borsani & Saintonge (2018) for neutral gas 
outflows of similar galaxies. 

Additionally, we also investigate the evolution of the out- 
flow EW, (geometry-independent) mass loss rate and mass load- 


ing factor (shown by the colour-coding of detections in the first 
three panels of Figure 7). Unlike in Section 4.2, the adopted N(H) 
value here is the mean column density per pixel. We find the 
first two quantities correlate positively with the outflow's posi- 
tion along the resolved MS: spaxels higher up the MS produce 
stronger outflows with significantly larger EWs and mass loss 
rates, consistent with the notion that more intense star formation 
activity drives stronger outflows. The mass outflow rates range 
from -3<log Mou/Moyr ! <0 with Hi column densities 18.6<log 
N(H)/em ? «22.0 and median values of log Mout ~-1.0 Moyr! 
and log N(H)~20.32 cm ?, respectively. However, the picture is 
not as clear cut for the mass loading factor, which does not appear 
to display significant evolution as a function of MS position (al- 
though there is tentative indication of higher factors further up the 
MS). The values here range from -0.4<log 7 <1.5, with a median 
mass loading factor of log 7 ~0.6. 

From the right panel of Figure 7, we find that outflows detec- 
tions are found predominantly in absorption and the lower limit of 
their detections corresponds to regions of reduced ISM absorption 
(or P-Cygni profiles), where the Nap profile transitions into emis- 
sion. Profiles of pure emission generally do not display any outflow 
signatures. The total Nap profile also shows a similar transition 
across the plane as observed in the SFR-M. plane with 3" fibers 
(Roberts-Borsani & Saintonge 2018; Concas et al. 2019): absorp- 
tion at high stellar mass and SFRs transits to emission at low values 
of the same quantities. However, the dichotomy of the two profile 
types with log Xsrg is more evident than with log X«, log SFR, or 
log M., and as such likely assigns Xgpg as the main regulator of 
the ISM profile. 

It is important to note that the evolution in outflow prevalence 
and properties seen over the local Xspg-3« plane also corresponds 
to the same radial evolutions seen in Figure 6, given that the evolu- 
tion of these particular galaxy properties also evolve with galaxy 
radius (i.e., both Xspg and X« decrease with increasing radius). 
Thus, the strongest outflow detections are found predominantly in 
the central regions of galaxies. 


6 OUTFLOW CORRELATIONS WITH GALAXY 
PROPERTIES 


As hinted at by Figure 6, the evolution of outflow properties is 
likely tied to one or several underlying galaxy properties. This has 
often been investigated for integrated galaxy quantities (such as to- 
tal SFR, stellar mass and Xgpg; Chen et al. 2010), however such 
investigations generally rely on single or stacked spectra of the cen- 
tral regions of a galaxy, making the isolation of the regions of in- 
terest challenging. 

However, due to the power of IFU spectroscopy, we are in the 
fortunate position of being able to identify and separate individ- 
ual spaxels (and therefore kpc-scale regions) in individual galaxies 
corresponding to a specific range of a given global galaxy prop- 
erty, thereby removing many of the intrinsic correlations that ex- 
ist between a given property and e.g., galaxy radius. We therefore 
perform this analysis for key galaxy quantities such as Xgppg, Ux, 
specific SFR (sSFR), Ay, and D(4000), which for each galaxy we 
divide into galactocentric radial bins in order to eliminate under- 
lying correlations between properties, and stack over all outflow 
galaxies. We present this analysis in Figure 8. Our choice of galaxy 
properties is motivated in part due to the availability of the tracers 
from optical spectra, but more importantly due to their inferred in- 
fluence on outflows from integrated studies: star formation-related 


MNRAS 000, 1—20 (2019) 


Resolved Outflows in Star-forming Galaxies 9 


OF @ outflow 


XL 
asg ce 
X LJ 


Ex 


‘outflow EW [A] 


[— HEN 


0.0 


log n [Mout/SFR] 


log Mo, [Moyr™*] 


log Zser [Moyr-!kpc-?] 


—3.00 -1.75 —0.50| 
10 7 
log Z- [Mckpc?] 


—0.50 0.25 


Figure 7. The resolved Xspg-2« plane of star-forming MaNGA galaxies with stacked bin limits (gray lines) and outflow detections (diamonds) overplotted. 
Non-detections are marked by a cross and detections are colour-coded according to the mean outflow EW (left panel), neutral gas mass outflow rate (middle 
left panel), and local mass loading factor (middle right panel). Outflows are prevalent in regions of high star formation activity and stellar mass surface density, 
with an EW and mass outflow rate that follows a similar pattern. For the mass outflow rate and loading factor, no outflow geometry is assumed. The right panel 
is the same as the other panels except colour-coded by the total Na D EW, as measured from the stacked spectrum. The dashed navy line indicates the resolved 
MS relation for H n regions of star-forming galaxies from Hsieh et al. (2017). 


quantities are invoked as the main drivers for outflows in the ab- 
sence of an AGN (e.g., Heckman et al. 2000; Veilleux et al. 2005; 
Chen et al. 2010), whilst both stellar mass and dust are influential 
in regulating the escaping potential of outflows and the survivabil- 
ity of Nap respectively. The D(4000) break traces the age of the 
underlying stellar populations, and as such can provide a first or- 
der indication as to which stellar populations may be driving the 
outflowing gas. 

We begin by noting that outflows are detected over a 
large range of each galaxy property. The ranges span -2.25<log 
Espgp/Moyr !kpc ? S-0.25, -11<log sSFR/yrS-9, 7.5<Slog 
X«Mokpc ? <9.5, I.2X D(4000) 22.0, and 0SAy/mag<3, con- 
sistent with ranges presented in all of our stacked analyses so far. 
Most of the detections are found in the inner radii, with the detec- 
tion rate falling rapidly as a function of radius. 
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Figure 8. The correlations between outflow properties (Cf, N(H), Av, Mout and 7) as a function of galaxy properties as shown by stacks of individual spaxels within a given property range. The mass outflow rates 
and mass loading factors here are independent of outflow geometry. The main trends occur with star-forming quantities, although the strongest correlation is seen with Xsrg. 
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Next, we inspect the local covering fraction of the gas, Cy, 
which is an estimate of its clumpiness; low values indicate that 
the gas is clumpy and therefore not covering much of the back- 
ground continuum source, whilst values closer to unity indicate a 
more diffuse component. For SDSS-selected galaxies, small values 
(Cf «0.5) for Nap gas have been found, and this is attributed to the 
low ionisation potential (5.1 eV) of the atom, which requires sig- 
nificant self- and dust-shielding to survive. As such, a more com- 
pact and clumpy nature ensures it can survive harsher environments 
(e.g., from shocks and hard stellar radiation fields). We find signif- 
icant scatter in the covering fraction of the gas spanning the full 
range of allowed values, with no clear evolution as a function of 
galaxy property or galaxy radius. Given that these stacks are con- 
structed over a variety of radii and the estimates of Cy are velocity- 
independent measurements, the scatter is perhaps not surprising. 

We also look at the column density of the outflowing hydro- 
gen gas, since it is interesting to determine where the densest re- 
gions of the outflows reside. We find clear, positive correlations of 
the column density with log Xspg and log sSFR, and a minor corre- 
lation found with log X«. No correlation is found with the D(4000) 
break and only a tentative trend with Ay may be present in the in- 
ner radii of galaxies. The first three quantities all strongly correlate 
with the amount of cold gas present and the star formation activity 
of galaxy regions, therefore such trends are perhaps not surprising: 
higher concentrations of gas induce more star formation which pro- 
duce stronger outflows which expel more gas. This is particularly 
evident from the strongest correlation with log Xspgg which spans 4 
dex in column density and suggests the densest parts of the outflow- 
ing gas correlate directly with the densest regions of star formation 
and cold gas. 

In measuring the blueshift of absorbed gas, we report abso- 
lute outflow velocities of «50-200 km s^ , similar to those reported 
for normal galaxies at z «O0 (Chen et al. 2010; Sugahara et al. 
2017; Roberts-Borsani & Saintonge 2018; Rodríguez del Pino et al. 
2019), however we observe essentially no correlations of the out- 
flow velocity with any of the chosen galaxy properties. The evolu- 
tion of the outflow velocity with galaxy properties and redshift has 
long been debated, with some studies claiming a velocity evolution 
over SFR (Sugahara et al. 2017) and others showing little to no 
evolution (Martin 2005; Chen et al. 2010; Roberts-Borsani & Sain- 
tonge 2018; Rodríguez del Pino et al. 2019). Such debate is also 
subject to the manner of measuring outflow velocities: for instance, 
in stacked spectra at z ~0, Sugahara et al. (2017) find convincing 
evidence for a velocity evolution across SFR when measured at the 
maximal blueshifted velocity of the absorption, however this evo- 
lution disappears when using the central velocity of the outflow, 
as shown here. Evolution in outflow velocity is also found to be 
more prominent at higher redshifts, with larger values of SFR and 
XspR found to correlate with enhanced outflow velocities traced by 
neutral and ionised gas at z ~2 (Sugahara et al. 2017; Davies et al. 
2019). As such, our results are in agreement with outflow velocities 
found in normal galaxies at z ~0, but in contrast with what is found 
by studies at z ^2. 

Perhaps the most important quantities to look at, however, are 
the evolution of the (geometry-independent) mass outflow rate and 
mass loading factor as a function of galaxy property. Here, we are 
able to directly link the kpc-scale galaxy properties to the outflow 
by selecting spaxels with relevant galaxy properties. Using Equa- 
tion | and the mean column density per spaxel, we find log Mout 
correlates most strongly (based on the slope and range of the ob- 
served correlations) with quantities most closely related to star for- 
mation activity, however the strongest correlation is associated with 
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log Xgpg: the mass outflow rate increases rapidly from log Mout ~- 
4 at log Xggg «-2.25 to log Mout ~0 at log Xgpg ~0.25. Milder 
but significant correlations are seen in the inner regions of galaxies 
with log sSFR, log X+, and Ay, likely due to their own correlations 
with log Xsrg and the survivability of Nap, although these all have 
both shallower slopes and probe smaller ranges. 


Finally, we also show that the mass loading factor of the out- 
flows again correlates most strongly with star formation activity, 
with values ranging -1<log 7 <1 in the central regions, and gener- 
ally following the evolution of the mass outflow rate. These values 
and trends are similar to those found for outflowing ionised gas in 
star formation-driven outflows by (Rodríguez del Pino et al. 2019) 
with MaNGA DR2, who also report loading factors of 7 <1. We 
are also in partial agreement with reports from Davies et al. (2019), 
who determine strong evolution between 7 and Xspg from ionised 
gas at z ~2: we observe a virtually identical slope and report similar 
values of the loading factor, although the overlap in data is small. 
Thus, from this analysis we can can infer that (1) the mass outflow 
rate is driven primarily by the evolution of log Xspg and (ii) the 
loading factor is most closely tied to the evolution of the mass out- 
flow rate (and hence also log XgpR), and (iii) that in the absence 
of extreme AGN, outflows are predominantly driven by the star 
formation activity in galaxies, traced most strongly by its surface 
density. Whilst such trends are clear, as mentioned in Section 4.2, 
our ability to detect blueshifted Na p in absorption is highly depen- 
dent on the strength of the underlying continuum - and thus Xy, 
- in order to give rise to net Nap absorption profiles. As such, we 
note that some degeneracy may exist between the outflow detection 
limit for some galaxy properties and the associated strength of the 
continuum that drives the EW of the Nap profile. This consider- 
ation should be kept in mind when analysing the results - and in 
particular, the threshold detection values - from Figure 7 and Fig- 
ure 8, where some underlying bias towards higher values of Xy, 
may be present. 


In addition to the above, we also look at the evolution of both 
the Nap ISM EW and the blueshifted absorption EW as a func- 
tion of galaxy properties, shown in Figure 9. For the Nap ISM 
EW, strong positive trends become immediately clear, with the 
strongest trends again occurring with quantites related to star for- 
mation (Zsrg. sSFR, X). A tentative inverse relation is found with 
the D(4000) index (with older stellar populations pushing the Nap 
EW to negative values), but curiously no significant trend is found 
with dust, despite it being thought to regulate (in part) the surviv- 
ability of Nap in harsh environments. Once again, the strongest of 
the trends is found with Xgpg and we interpret the increase EW 
and Na p absorption as being due to the increased presence of cold 
gas with higher levels of Xspgg. We also note that the majority of 
Nap is found in absorption, particularly at higher values of Xspg. 
sSFR, X. and Ay, but transforms into net emission at lower val- 
ues. Interestingly, the transition of Nap from absorption to emis- 
sion in stacks of Xspg occurs right around our lower detection limit 
at 0.01 Mayr" !kpc"? and consistent with what we find in the 
right panel of Figure 7. The evolution of the outflow EW as traced 
by blueshifted absorption is similar but less pronounced than what 
we observe for the entire ISM line, in the inner regions of galaxies: 
positive trends are found with Xspg, sSFR and dust, with higher 
values corresponding to stronger outflowing gas, whilst log X. and 
X., D(4000) show no clear trends. Such findings are consistent with 
the conclusions derived above and support the notion that quanti- 
ties related to star formation are the main regulators of both the 
Nap EW and outflow EW. 
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Figure 9. The evolution of the Na p ISM EW (top row) and blueshifted absorption EW (bottom row) as a function of Xspg, sSFR, £., D(4000) and Ay. The 


colour scheme follows that of Figure 8. 


7 THEIMPACT OF OUTFLOWS ON HI GAS 
RESERVOIRS 


7.1 HlObservations 


One of the most important questions in outflow studies is whether 
they have some impact on the cold atomic and molecular gas reser- 
voirs of their host galaxies (i.e., can they quench the galaxy?). So 
far, we have inferred this from derivations of a neutral gas mass 
loading factor. However, an alternate approach is to examine the ef- 
fect of outflows on the integrated H1 contents of their hosts. To do 
this, we make use of H 121 cm single-dish observations of galaxies 
in our outflow+non-detection samples with the Green Bank (GBT) 
and Arecibo Telescopes. The GBT observations form part of the 
HI-MaNGA (Masters et al. 2019) program aimed at carrying out 
follow up Hı observations for MaNGA galaxies and the DR1 con- 
tains 331 galaxy observations to date. Of these, 33 are matched 
to our selected sample, however we also add an extra 53 galax- 
ies which were observed after the data release (Karen Masters and 
Dave Stark, private communication), bringing the MaNGA GBT 
observations used here up to 86 galaxies. Additionally, some of our 
sample also overlaps with observations from the blind ALFALFA 
survey (Haynes et al. 2018) conducted with the Arecibo telescope. 
55 galaxies from our sample are matched to ALFALFA and as such 
we are able to compile H r observations for a total of 141 galaxies in 
our combined outflow and non-detection samples, of which 34 fall 
under the outflow category and 107 in the non-detection category. 
We refer to this as our MaNGA-H 1 sample. 


7.2 Removing the Effects of Confusion and Baseline Issues 


Spectroscopic confusion in radio observations from single-dish fa- 
cilities is an important concern, and in particular for stacking and 
outflow studies, since real signal from nearby galaxies at similar 
velocities can mimic the signatures of outflowing gas. 

Although the rate of spectroscopic confusion is generally 
small (less than ~15% for the ALFALFA a.40 data release; Jones et 
al. 2015), given the large beam sizes of the facilities used to obtain 
the Hı data in this study (FWHM*9’ and FWHM3.5' at 21cm 
for the GBT and Arecibo, respectively) and the potential for false- 
positive outflow detections, this remains a particularly important 
consideration. 


Therefore, to ensure none of our spectra suffer from con- 
fusion, we use the MPA-JHU catalog to identify objects within 
3xFWHM of the beam used to observe each of our 141 Hı galax- 
ies (but further than 10" away from the target galaxy). Combining 
this with a velocity cut of +500 km s^! from the target source (rel- 
ative to the H velocity of the target if detected, or optical redshift 
in the case of non-detections), we flag 78/141 of our objects as 
containing one or more additional optical galaxies over the search 
area and potentially subject to confusion. The criteria used here are 
extremely conservative, however we believe a rigorous approach 
to avoid false-positive outflow detections is crucial, and given that 
~80% of the contaminating sources fall outside the FWHM of the 
beam (where the sensitivity drops from 50% to effectively 0%) 
combining the analysis with the beam sensitivity (see below) is 
likely to decrease this number significantly. 


We next estimate the Hı gas masses of the contaminants, us- 
ing gas fraction scaling relations (Catinella et al. 2018), and in par- 
ticular the relation between log My ;/M.. and log sSFR. However, 
given that the sensitivity of the beam drops dramatically past its 
FWHM and we probe an area a factor of 3 larger (effectively span- 
ning the full sensitivity of the beam) we multiply the derived H1 
masses by the beam sensitivity at their separation in order to obtain 
an effective Mg; and compare this to the target's H1 gas mass. If a 
given contaminant contributes less than 10% to the target’s H 1 mass 
(in the case of a non-detection we assume upper limits as calculated 
in Masters et al. 2019), it is considered to have a negligible effect 
on the spectrum (Jones et al. 2015). We flag a galaxy as confused 
if it has one or more surrounding contaminants (within the velocity 
cut defined above) that contribute >10% of the target's measured 
Hı gas mass. 32/141 of our MaNGA-HI sample are flagged as con- 
fused and discarded from the analysis. 


A second important consideration in our stacking analysis is 
the possible effects resulting from badly removed baselines, which 
can result in artificially low or high fluxes over various regions of 
the spectrum. To ensure our spectra are free from such effects, we 
visually inspect each of the remaining 109 galaxies that are free 
from confusion and determine whether a baseline correction is re- 
quired. 6 galaxies are found to have major baseline issues, with 
an additional 7 found with minor issues. For these galaxies, we fit 
the baseline-unsubtracted spectrum with a polynomial (generally 
of 3rd or 4th order) to the baseline and subtract this from the spec- 
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trum. During this procedure we consider only the regions within 
+1000 km s^! of the target galaxy velocity, since fitting an accu- 
rate baseline across the entire spectrum can be challenging and our 
region of interest is largely confined to those velocity limits. Only 
2 galaxies have major baseline issues which we are unable to cor- 
rect for and an additional target contains obvious spectral artefacts, 
and are therefore discarded from the analysis, leaving a total of 106 
galaxies free from confusion and baseline issues, which we refer to 
as our MaNGA-HIcorr sample. 


7.3 Control Sample and Stacking Procedure 


In this analysis, we wish to determine whether there is significant 
difference between the H1 gas reservoirs of galaxies with Nap out- 
flows and those without. Thus, for each of the outflow galaxies in 
our MaNGA-HIcorr sample, we select a control galaxy without 
Nap outflows which is matched in its position on the SFR-M.. plane 
and in inclination. We allow a difference of +0.2 dex in log SFR 
and log M,, as well as 20 degrees in inclination. Of our MaNGA- 
HIcorr sample, we can successfully identify control galaxies for 
17 of our outflow galaxies. 

In order to construct mean H i spectra, we stack in “gas frac- 
tion units", that is to say we multiply each spectrum by the stan- 
dard conversion factors necessary to convert to an H 1 gas mass, de- 
scribed in Equation 6 and divide by the stellar mass of the galaxy. 


2 
Muy /Mo = 2.356 x (S (=) (6) 
Mpc Jy 

All scaled spectra are first shifted to the rest frame prior to 
being normalised by the FWHM of the spectrum, then interpolated 
over a common velocity grid, and finally added to the stack irre- 
spective of whether they display a clear Hı detection or not. The 
stack is then averaged to produce a mean gas fraction spectrum. 
We further adopt a Monte Carlo approach and repeat this process 
100 times, each time with a random sample of 14 outflow galax- 
ies (~80% of our outflow H1 sample) and a different control sam- 
ple, to ensure our results are not biased by a particular selection of 
galaxies. The final spectra are then taken as the mean over the 100 
iterations and shown in Figure 10. 

The resulting velocity-normalised spectra display significant 
differences between themselves, with outflow galaxies displaying 
reduced fluxes compared to their non-detection counterparts, sug- 
gestive of non-negligible variation in the Hı gas mass fractions. 
The mean FWHM is 209 km s^! for the outflow galaxies and 205 
km s^! for the non-detection galaxies, with integrated Hı gas frac- 
tions of log fur--0.90620.001 and log fur--0.75720.001, respec- 
tively. We also present an inset plot in Figure 10 displaying the dis- 
tributions of gas fractions of the individual galaxies going into the 
stack, over all Monte Carlo iterations. The two distributions are dis- 
tinctly offest, with means of log far=-0.91 and log fur--0.76 (out- 
flow detections and non-detections, respectively) consistent with 
the integrated spectra, and spreads characterised by standard devi- 
ations of log o, ~-1.79 (outflow detections) and log ag, ~-1.47 
(non-detections). The two distributions reinforce the apparent lack 
of similarity between the two stacked spectra and the combination 
of the two is suggestive of a non-negligible difference between the 
Hı gas fractions of the outflow-selected galaxies and those without 
outflows. 

Although causality is challenging to establish, we have been 
careful to match the outflow-detection and non-detection galaxies 
in several key properties that could influence the observed H1 gas 
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Figure 10. The mean velocity normalised, gas fraction spectra of galax- 
ies with signatures of optical outflows (blue) and those without (gray) and 
their reported gas fractions and median FWHM. The stacks are created via 
a Monte Carlo approach of 14 outflow galaxies and a control sample of non- 
detection galaxies. The presence of Na p -selected outflows does not appear 
to significantly influence the H I gas reservoirs of their host galaxies com- 
pared to their non-detection counterparts. The inset plot shows a histogram 
of the average gas fraction of a given iteration going into the final stack. 


fraction, to ensure the primary difference is the detection of an op- 
tical outflow. Thus, the different shapes of the profiles (a double 
horn profile for outflow galaxies and a Gaussian-like profile for the 
non-detection galaxies) could hint at blowing out of H1 gas at the 
centres of the galaxies by the outflows (i.e., at normalised velocities 
of ^0 km s-l), where the outflows are strongest (see Section 4.2). 
Furthermore, the higher Hı gas fractions at virtually all velocities in 
the non-detection stack are suggestive of some gas replenishment 
and could provide indirect evidence for the interplay between out- 
flows and inflows - i.e., a galactic fountain scenario. Whilst, such 
a scenario is a distinct possibility, we also note that such features 
could also reflect some transition of neutral Hı gas into molecular 
or ionised gas, in which case more molecular gas may be available 
for star formation in the galaxy centres and could lead to the forma- 
tion of the winds. However, distinguishing between and confirming 
such scenarios is beyond the capabilities of the current data sets. 


8 DISCUSSION 
8.1 Strong Correlations With Xspg 


Throughout each of the previous Sections, we have seen that Xspg 
plays a key role in regulating the output of outflows in normal 
galaxies at z ^0. Whilst other quantities (e.g., sSFR, X. and Ay) 
also correlate with outflow quantities, these appear secondary to the 
dominant correlations found with Xgrg and likely exist due to their 
own correlations with Xspg. Similar correlations have also been 
observed in other studies, both at low and high redshifts (Chen et 
al. 2010; Newman et al. 2012; Roberts-Borsani & Saintonge 2018; 
Davies et al. 2019). For galaxies at z ~2, Newman et al. (2012) 
suggested the possibility of a critical “blow out" Xspg, where star 
formation feedback is able to generate enough pressure to per- 
pendicularly break out of the dense gas of the disk in the form 
of a momentum-driven outflow. By assuming a baryon dominated 
galaxy disk that sits in pressure equilibrium, one can test this via 
simple assumptions that equate the weight of the disk gas to the 
pressure exerted outward by star formation feedback, and deter- 
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mine a minimum Xgpg threshold above which pressure from the 
feedback exceeds the weight of the gas. 

Following Equations 1 and 7 from Ostriker, & Shetty (2011), 
the weight of the disk can be expressed as 


(7) 


where Yeas is the cold gas density and G is the gravitational 
constant. The vertical momentum flux injected into the ISM by stel- 
lar feedback can be described as 


P= 2 Yepp, (8) 
Mx 

where (p../m..) is the mean radial momentum injected into the 
ISM per unit mass of stars formed. As described by Equation 2 of 
Newman et al. (2012), one can rearrange these two equations and 
insert a dependence on gas fraction, fg and stellar mass density to 
derive a minimum Xsrg threshold. This is particularly useful given 
that such quantities are well constrained for local, normal galaxies 
and by the data presented here. As such, our final equation is 


nG fg o 
XSER,thresh = p,m) (9) 


Assuming fg ~0.5, (p+/m)~1000 km sland Zgas ~500-1000 
Mopc ? for their sample of normal SFGs at z ~2, Newman et 
al. (2012) found a critical threshold of ^1 Moyr"!kpc ?, in good 
agreement with their observations of ionised outflows. Assuming 
typical values (fg ~0.07, Esas ~500) for z ~0 normal SFGs, this 
threshold is lowered down to ~0.1 Mayr" !kpc"?, in agreement 
with the canonical value typically assumed as the minimum thresh- 
old to launch outflows (Heckman 2002; Veilleux et al. 2005; Os- 
triker, & Shetty 2011). 

However, this threshold is typically observed in more turbu- 
lent, starburst systems and not in normal SFGs along the MS. 
Here, we expand on this by using values directly inferred from 
our MaNGA data set. Assuming a cold (molecular) gas fraction for 
high mass galaxies of fọ ~0.04 (Saintonge et al. 2017) and a mean, 
galaxy-wide mass surface density of log XJMokpc ? ~8.24 over 
all of our sample of 376 outflow+non-detection galaxies, we de- 
rive a critical threshold of Zgpg ~0.01 Moyr-!kpc^?, in agreement 
with the results presented here and in Roberts-Borsani & Saintonge 
(2018). Additionally, we also note that the mean galaxy-wide Xsrg 
for our outflow selected galaxies is 0.02 Mo yr"! kpe~? whilst that 
of the non-detection galaxies is ~0.01 Moyr~'kpe~?, sitting above 
and on our derived threshold, respectively. 

For reference, we compare this value to the average SFR sur- 
face density of the Milky Way, which does not harbour a strong out- 
flow (Fox et al. 2019). Assuming a SFR of 1.65 Mo yr! (Bland- 
Hawthorn & Gerhard 2016, and references therein) and a disk ra- 
dius of 10 kpc (Olsen et al. 2015, and references therein), the area of 
the disk is 314 kpc? and thus the SFR surface density is calculated 
to be spr ~0.005 Moyr. 4kpc ?, below our derived threshold. 
Such values can provide a useful comparative framework, however 
similar analyses over the central regions of the Milky Way (e.g., 
the Central Molecular Zone or the Galactic Centre) could provide a 
more informative comparison given they likely provide the regions 
with the strongest outflowing gas, whilst the star formation history 
of these regions should also be taken into account. 


8.2 Star Formation Histories of Outflow Hosts 


Given that star formation quantities appear to be the primary drivers 
of our selected outflows, it is instructive to look at the star for- 
mation histories associated with the host galaxies which give rise 
to the outflows. As such, we use the D,,(4000) and Hô; indices 
provided by the MaNGA DAP, whose combination can reveal and 
distinguish between continuous and bursty star formation histo- 
ries (Kauffmann et al. 2003). The D,(4000) index measures the 
strength of the 4000 À absorption break, which is a discontinuity 
in the optical spectrum of galaxies due to a variety of absorption 
features primarily arising from ionised metals in the atmospheres 
of stars. The index is linked to the age of the galaxy's stellar pop- 
ulations, given that young, hot stars ionise the surrounding metals, 
leading to a decreased opacity, whilst this is not the case for old, 
metal rich stars which display considerably more metal absorption. 
As such, a large 4000 À break is indicative of older stellar popu- 
lations, whilst a smaller one is indicative of younger stellar popu- 
lations and a generally increasing break is found with time, as the 
stellar populations evolve. The index itself is measured as the ratio 
of the average continuum flux in regions blue and redward of the 
break, with several sources in the literature adopting differently- 
sized regions (broad or narrow). The advantage of adopting a more 
narrow-sized region either side of the break is that the ratio is less 
sensitive to the effects of dust. As such, we use the narrow-band 
definition adopted by MaNGA (3850-3950 and 4000-4100 À for 
the blue and red continuum bandpasses, respectively; Westfall et 
al. 2019). 

Hó absorption at 4101 A, on the other hand, arises from 
Balmer absorption in the atmospheres of stars and is an indication 
of the timing of a recent burst of star formation. Beginning with 
weak intrinsic absorption by hot O and B stars (characterised by 
short lifetimes of a few 10 Myrs) at the time of the burst, the ab- 
sorption increases monotonically over time until ~300 Myr where 
it reaches its peak. The rapid increase and peak of the evolution is 
due to the deaths of the OB stars and domination of the optical light 
by late-B to early-F stars (characterised by longer lifetimes of <1 
Gyr). After the peak absorption at ~300 Myr, the Hó absorption 
rapidly decreases as the A and F stars die off. As such, Hó absorp- 
tion is a measure of recent bursts of star formation that ended only 
~0.1-1 Gyr (Kauffmann et al. 2003), and its index (HóA) is mea- 
sured using a central bandpass marked by two pseudo-continuum 
bandpasses either side of the line (4083-4122, 4041-4079, 4128- 
4161 À for the main, blue and red bandpasses, respectively; West- 
fall et al. 2019). 

The two indices are largely independent of reddenning effects, 
however metallicity can play a role in regulating them at later times. 
However, as shown by Kauffmann et al. (2003), metallicity only be- 
comes an important consideration for both tracers at ages of »10? 
years after a burst (corresponding to rough values of D,4000>1.5 
and HóA «3-4). 

To conduct our analysis, we use the average of the two tracers 
within a galactocentric radius of 1.5" for each of our outflow and 
non-detection galaxies and compare these to results in the MPA- 
JHU catalog as measured by the SDSS DR7 3’’-diameter fiber. The 
SDSS relations are taken by binning star-forming galaxies selected 
by the same MaNGA selection criteria described in Section 3 in 
steps of AD,,(4000)=0.2 and taking the mean and median values, 
before fitting a 5th order polynomial. We present these fits and data 
points in Figure 11. 

We find our data points span a wide range over both indices 
and generally sit slightly above the mean (and median) SDSS rela- 
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Figure 11. The star formation histories as traced by the 3"-diameter stacked 
Ho, and D, (4000) indices of our MaNGA DR15 sample of galaxies shown 
to host Na p outflows (blue points) and those that do not (gray points). The 
error bars represent the standard deviation of values going into the stack for 
each galaxy. These points are compared to the mean (dashed orange line) 
and median (dashed red line) of SDSS values from the MPA-JHU catalogs. 
The inset plots show the distributions of AHóA (Hó A, MaNGA-HÓ A,spss: 
top) and HóA (bottom) for outflows (blue lines) and non-detections (gray 
lines), as well as their cumulative distribution functions and associated two- 
sample Kolmogorov-Smirnov statistic. 


tion, regardless of whether they are selected to have Nap outflows 
or not. For the outflow galaxies, ~63% of the data points sit above 
the mean SDSS relation and ~37% below, suggesting the majority 
of the data points have experience higher bursts of star formation 
in their recent past. However, comparing this to the non-detection 
data points, which reveal a slightly higher percentage of ~70% sit- 
ting above the SDSS relation and ~30% sitting below, we find very 
similar values. This is in contrast to IFU observations of ionised 
tracers in edge-on galaxies found by Ho et al. (2016), who ob- 
serve 8046 of their outflow galaxies and 5046 of their non-detection 
galaxies to sit above the median SDSS relation, albeit with sam- 
ple sizes a factor of ~5 and ~12 smaller than those used here, re- 
spectively. We expand on this by comparing the distribution and 
mean values of our data points with the corresponding SDSS val- 
ues (AHóA-HÓó A MaNGA-HÓ A, spss), inferred from our fit relation. 
We present this in the top inset plot in Figure 11. The two distribu- 
tions are very similar, with outflows characterised by a mean dif- 
ference of 0.27 and standard deviation of 0.78, and non-detections 
with a mean of 0.32 and identical standard deviation. This is rein- 
forced by near identical cumulative distribution functions (CDFs) 
and a low Kolmogorov-Smirnov statistic of 0.11, suggesting little 
difference between the two distributions and supporting the notion 
that both the outflow and non-detection samples have both expe- 
rienced more recent star formation bursts compared to their SDSS 
counterparts. Such similarity extends to the two MaNGA samples, 
which we compare via their HóA distributions in the bottom inset 
plot in Figure 11. The two distributions are once again quite sim- 
ilar, with means of 4.02 and 3.56 and standard deviations of 1.42 
and 1.82 for the outflow and non-detection samples, respectively. 
The two CDFs and a K-S statistic of 0.17 reinforce this similar- 
ity, albeit with a slight offset for outflow galaxies towards higher 
Ho, values. We do note, however, that the distribution of the non- 
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detection galaxies does fall down to HóA «3-4, where the tracer 
becomes uncertain due to metallicity effects, and the differences in 
mean values is attributed to this extension. 

The similarities between the two MaNGA samples suggest 
both the outflow-detected and non-detected samples may be under- 
going similar episodes of short bursts of more intense star forma- 
tion that may be capable of launching outflows. If this is the case, 
both samples would display similar values of HóA and D;,(4000), 
although one may also expect more of the non-detection galaxies to 
display signatures of optical outflows. Thus, the lack of (detected) 
outflow signatures in the non-detection sample is perhaps surpris- 
ing. Following on from this, whilst our results are in agreement 
with the ionised gas results of Ho et al. (2016) in that the majority 
of our outflow data points lie above the mean SDSS relation, the 
considerable similarity between our larger distributions of outflow 
and non-detection samples is in contrast to theirs and suggests that 
(i) our two MaNGA samples may be undergoing similar episodes 
of star formation and (ii) that the timing of the last burst of star 
formation may not be a crucial parameter in driving (weak) Nap 
outflows. 


9 SUMMARY AND CONCLUSIONS 


We use the SDSS-IV/MaNGA DR15 data set and spectral stacking 
techniques to constrain the main kpc-scale properties that give rise 
to and regulate neutral outflows in star-forming galaxies. We use 
a sample of 405 z ~0 face-on galaxies along the galaxy MS to 
determine the detection fraction, galactocentric radial profile and 
kpc-scale properties of galaxies with signatures of blueshifted Nap 
absorption. Our main findings can be summarised as follows: 


e Out of 390 useable galaxies in our sample, the stacking of 
the central 0.25 Re regions of the galaxies reveals 78 objects 
with signatures of outflowing gas, 14 objects with signatures of 
inflowing gas, and 298 galaxies with no significant detection of 
blue/redshifted gas. Galaxies with signatures of outflows and 
inflows show considerably higher values of X.. and XspR compared 
to their non-detection counterparts. 


e Stacking as a function of deprojected galactocentric radius, 
we find detections of outflowing gas out to 1 Re and compare these 
to detections as observed by an 3’’-diameter fiber. Derivations 
of a mass outflow rate show a rapidly declining trend with 
galactocentric radius, with a range -2.1<log Mout <-0.9 and best 
fit slope of -1.2. A near identical trend is found with the mass 
loading factor, which also decreases as a function of galaxy radius, 
demonstrating that outflows are at their strongest in the central 
regions («0.5 Re) of galaxies. The range of reported mass loading 
factors is -1.8<log 7 -0.4, consistent with previously reported 
values in relatively normal galaxies at z ~0, and a slope of -1.4. We 
provide prescriptions for both of these trends for use in simulations. 


e Signatures of outflowing gas are found along and above the 
resolved star-forming MS, in parameter space above log X. 27.5 
Mokpe~? and log Uspp Z-2 Moyr !kpc 2, similar to trends 
found along the galaxy MS. A comparison of the outflow EW, 
mass outflow rate and mass loading factor shows an increase in 
each of these values as one moves up the resolved MS, suggestive 
of stronger outflows in regions of increased star formation activity. 


e By stacking 257,500 individual spaxels associated with a 
variety of galaxy properties (Xspr, X., specific SFR, D(4000) 


16  Roberts-Borsani et al. 


and Ay), we find significant positive correlations between out- 
flow properties (namely, gas column density, mass outflow rate 
and mass loading factor) and all of the aforementioned galaxy 
properties, except D(4000). The strongest correlations are found 
with log Xsgg and we extend the lower limit of detections down 
to log Xsgg -2 Moyr !kpc?, about an order of magnitude 
lower than canonical values. We suggest that this is due to a 
minimum Xsrg threshold of Xspg 40.01 Moyr-!kpc? necessary 
for star formation feedback to break out of the dense gas in the disk. 


e Using Hı follow up observations of a sample of 106 MaNGA 
galaxies, we compare the cold Hı gas reservoirs of galaxies se- 
lected to have Nap outflows to control samples of galaxies without 
Nap outflows, matched by their position in the SFR-M, plane 
and inclination. Stacking the Hı spectra, we find considerable 
differences between gas reservoirs of outflow and non-detection 
galaxies, suggesting the presence of optically-selected outflowing 
gas could potentially have some (minor) effect on the cold gas 
reservoirs, although combining this with our Nap mass loading 
factors strongly suggests that it is not enough for normal galaxies 
to be quenched by ejective feedback. 


e Finally, we compare the star formation histories of our outflow 
and non-detection galaxies, using the central (3’’) stacked HóA and 
D,,(4000) indices. In comparing to average SDSS DR7 values, we 
find our sample of MaNGA galaxies have slightly elevated values 
of HóA compared to the SDSS galaxies for a given D,,(4000), how- 
ever, crucially we find virtually no difference between the values of 
our outflow sample and non-detection sample, suggesting that the 
timing since the last burst of star formation is not a crucial param- 
eter in driving Nap outflows. 


The arrival of large IFU surveys has greatly aided studies of 
local galaxy evolution and, in particular, studies of outflows. The 
kinematics and key properties of outflows and their hosts are be- 
ing studied in unprecedented detail and the small, kpc-scale pro- 
cesses that give rise to the large scale processes are being con- 
strained. Here, we have used samples of star-forming galaxies and 
>275,000 galaxy resolved spectra to determine the prevalence of 
outflows on kpc-scales, their power and quenching potential across 
a variety of galaxy regions and properties and determined whether 
any significant impact on the H1 gas reservoirs in seen. However, 
still missing are crucial constraints on the multiphase nature of 
outflows in normal, MS galaxies at z ~O in order to infer their 
true quenching potential. Verifying this will require rest-frame opti- 
cal, (sub)millimetre and radio spectra over statistically meaningful 
samples of galaxies with which to probe and compare the different 
outflowing gas phases, with the aim of deriving total mass loading 
factors over large regions of parameter space. Large surveys such 
as the SDSS, xCOLD GASS and ALFALFA provide unique oppor- 
tunities to do this, whilst the arrival of optical and radio surveys 
such as MaNGA, SAMI and PHANGS will allow similar analyses 
to be achieved on resolved scales. 
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APPENDIX A: GLOBAL PROPERTIES OF FLOW 
GALAXIES 


We present here the main global properties of MaNGA DRI5 
galaxies identified to host either an outflow or inflow in the central 
regions of the galaxy, whose properties were used in the sample se- 
lection. Quantities with superscript a mark quantities derived from 
the NSA catalog, whilst quantities with superscript b are derived 
from the Pipe3D catalog. SFRs and stellar masses are integrated 
quantities, whilst the line ratios refer to measurements in the cen- 
tral 2.5" of the galaxy. 
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Table Al: The global properties of galaxies identified to host outflows in their 
central regions. 


Plate ID  IFU RA? DEC? z? log M, log SFR®  log[Om]/Ho? log [NuJ/Hp? i? Type 
[deg] [deg] [Me] [Moy] [deg] 
7815 3702 317.90320 11.49694 0.02938 10.47 0.43 -0.73 -0.36 27.49 outflow 
8077 3702 41.84637 0.05876 0.02478 10.12 -0.10 -0.59 -0.37 23.92 outflow 
8085 3704 51.15552 -0.44402 0.03684 10.72 0.66 -0.74 -0.44 44.23 outflow 
8134 1901 113.40018 45.94338 0.07665 10.98 1.29 -0.63 -0.49 27.00 outflow 
8146 1901 117.05387 28.22509 0.02706 10.38 0.42 -0.69 -0.43 22.85 outflow 
8147 6103 119.03950 26.87578 0.06161 11.06 1.00 -0.72 -0.48 36.35 outflow 
8149 3704 120.20722 27.50036 0.01734 10.30 0.62 -0.55 -0.43 32.84 outflow 
8155 3703 54.090596 -0.88266 0.02358 10.18 0.27 -0.11 -0.52 25.001 outflow 
8253 3703 157.34358 43.17059 0.02726 10.32 0.39 -0.72 -0.47 45.99 outflow 
8254 1902 163.21911 43.42840 0.02424 10.09 0.02 -0.81 -0.37 27.02 outflow 
8262 9102 184.55357 44.17324 0.02453 10.66 0.72 -0.80 -0.34 39.23 outflow 
8309 6104 211.78062 52.96374 0.04333 10.21 0.33 -0.50 -0.53 38.35 outflow 
8311 3703 205.01217 23.14297 0.03168 10.84 0.60 -0.40 -0.34 34.63 outflow 
8312 12702 245.27087 39.91739 0.03203 10.70 0.52 -0.51 -0.30 43.71 outflow 
8312 12703 247.20906 39.83509 0.03585 10.83 0.73 -0.46 -0.36 16.04 outflow 
8313 12701 239.49180 41.79259 0.03547 10.78 0.72 -0.52 -0.45 47.27 outflow 
8313 12702 240.67742 41.19726 0.03334 10.73 1.06 -0.69 -0.49 48.87 outflow 
8317 6102 194.92503  43.75318 0.05781 11.13 1.16 -0.36 -0.39 35.55 outflow 
8329 6104 213.11048 45.69041 0.02700 10.54 0.73 -0.43 -0.32 31.33 outflow 
8332 3702 207.87281 43.80643 0.03331 10.56 0.67 -0.28 -0.39 39.89 outflow 
8333 3701 214.34452 42.28797 0.06861 10.84 0.72 -0.40 -0.34 27.82 outflow 
8446 6104 206.54141 36.50439 0.05524 10.54 0.84 -0.52 -0.46 37.46 outflow 
8454 12702 153.04902 45.14216 0.07643 11.13 1.40 -0.33 -0.36 44.77 outflow 
8454 6104 154.82544  45.81359 0.06348 11.07 1.12 -0.47 -0.33 34.15 outflow 
8455 3701 157.17946 39.83888 0.02929 10.34 0.60 -0.58 -0.38 41.65 outflow 
8465 12705 198.23632 47.45663 0.02814 10.81 0.11 -0.44 -0.34 47.01 outflow 
8465 9102 198.18916 46.93499 0.02768 10.40 0.12 -0.83 -0.47 33.48 outflow 
8486 6101 238.03955 46.31979 0.05887 10.88 0.60 -0.44 -0.33 43.53 outflow 
8550 12703 247.67443 40.52939 0.02981 10.49 0.05 -0.65 -0.33 31.70 outflow 
8550 9102 247.20906 39.83509 0.03585 10.81 0.73 -0.46 -0.35 16.04 outflow 
8551 12705 233.94092 44.83480 0.02964 10.59 0.79 -0.56 -0.40 26.90 outflow 
8554 3704 184.62322  35.62226 0.03490 10.44 0.35 -0.60 -0.45 29.22 outflow 
8588 6101 248.45676 39.26321 0.03176 10.88 0.67 -0.41 -0.34 16.10 outflow 
8600 1901 242.58522  43.00964 0.02521 10.19 0.48 -0.16 -0.58 46.99 outflow 
8611 3702 261.46394 60.19412 . 0.02905 10.51 0.13 -0.53 -0.35 41.79 outflow 
8616 3703 322.51023 0.46420 0.13506 11.15 1.21 -0.56 -0.38 28.05 outflow 
8618 3704 318.86229 9.75782 . 0.07023 10.94 0.82 -0.59 -0.41 10.25 outflow 
8624 9102 263.89258 59.88992 0.02836 10.55 0.27 -0.54 -0.36 48.19 outflow 
8626 12703 265.10878 56.70125 0.06928 10.92 0.66 -0.34 -0.33 41.12 outflow 
8626 3703 264.66254 56.82424 0.02942 10.33 0.51 -0.45 -0.44 38.56 outflow 
8715 9101 119.91012 51.79236 0.05338 10.72 0.63 -0.63 -0.35 48.18 outflow 
8725 3703 125.82096 45.89686 0.05373 10.64 0.59 -0.59 -0.34 44.09 outflow 
8725 6102 125.45904 45.51964 0.05417 10.66 0.87 -0.66 -0.33 22.84 outflow 
8940 1902 120.98713 25.48072 0.07278 11.10 1.29 -0.60 -0.38 29.88 outflow 
8945 3702 173.36190 47.28673 0.04555 10.05 0.65 0.07 -0.67 29.13 outflow 
8977 9102 119.20025 33.25448 0.06756 10.86 1.16 -0.47 -0.46 35.26 outflow 
8992 12704 172.63126 51.45594 0.02559 10.55 0.48 -0.77 -0.52 40.50 outflow 
9026 3703 249.34435 44.34712 0.03069 10.48 0.66 -0.87 -0.51 19.68 outflow 
9028 9101 243.02559 28.49887 0.05318 10.83 0.84 -0.84 -0.45 27.34 outflow 
9029 12702 247.25172 41.28426 0.03189 10.80 0.57 -0.66 -0.49 33.56 outflow 
9029 1902 246.48724  41.52207 0.04272 10.03 0.01 -0.74 -0.32 22.29 outflow 
9033 1902 223.21854 45.23480 0.13265 11.16 1.27 -0.55 -0.41 37.88 outflow 
9033 3703 222.93490 47.86576 0.04988 10.56 0.53 -0.45 -0.35 27.50 outflow 
9035 6102 236.65444 43.68365 0.10117 11.03 1.18 -0.55 -0.40 47.29 outflow 
9041 9102 235.94192 28.41521 0.03275 10.66 0.67 -0.55 -0.46 33.55 outflow 
9042 12703  235.15268  28.51244 0.03275 10.71 0.23 -0.62 -0.37 46.84 outflow 
9044 6101 230.68698 29.76960 0.02292 10.24 0.39 -0.78 -0.47 33.96 outflow 
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